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ABSTRACT 


Generalized analytic results are derived to give a complete description of an open- 
loop minimal control requirement to maintain stability for a rigid-body vehicle configura- 
tion. The results are easily applicable to both pitch and yaw control planes at any flight 
time of interest by computing some vehicle dependent parameters and evaluating the given 
equations. The results show that, with good control system design, the deflection re- 
quirements can be reduced to about 56 percent of trim deflection. At the same time, be- 
cause of the method of control, the maximum angle of attack and bending moments can 
also be substantially reduced. 
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STABILITY ANALYSIS AND MINIMUM THRUST VECTOR CONTROL 
REQUIREMENTS OF BOOSTER VEHICLES DURING 
ATMOSPHERIC FLIGHT 
by Janos Borsody and Fred Teren 
Lewis Research Center 

SUMMARY 

Open-loop stability boundary conditions are established for a large class of vehicle 
configurations. The results can be used for preliminary design purposes. The stability 
boundary is defined as the minimum thrust-vector deflection required to maintain a stable 
flight. This is accomplished through a fixed-time analysis, which assumes that the basic 
vehicle parameters remain constant in the time interval of interest. A rigid-body con- 
figuration is assumed with uncoupled pitch and yaw control planes. The equations are 
derived for the pitch plane, but can be applied to the yaw plane with minor modifications. 
The calculus of variations is used to derive an optimum thrust -vector deflection angle 
profile with the smallest maximum absolute value, and this control action is used, with 
minor modifications, to predict the stability boundaries presented. 

The results obtained show that the thrust-vector deflection requirements can be re- 
duced to 56 percent of trim requirement with a good closed-loop control system design. 
The analysis does not include the design of such a system. Attitude response is pre- 
sented for different wind disturbances and thrust-vector deflection profiles. These re- 
sults indicate rapid convergence and divergence of attitude error for control capability 
above or below the stability boundary value, respectively. 


INTRODUCTION 

Most present booster vehicles use thrust vectoring for control. This may be accom- 
plished in one of several ways, but the most common method is to gimbal the engine or 
engines. For liquid-fueled rockets, this method is used almost exclusively. However, 
thrust -vector deflection (TVD) of solid-fueled motors presents problems not encountered 


with liquid-fueled engines. It is not feasible to gimbal the entire motor, and, if only the 
nozzle portion is gimballed, severe seal and actuator problems are encountered. Sec- 
ondary fluid injection has been shown to be an effective steering technique, but, for many 
vehicle requirements, it leads to increased complexity, large injectant requirements, 
and reduced payload capability. Furthermore, any steering which is accomplished by 
generating side forces near the base of the vehicle introduces large bending moments in 
the vehicle structure in flying through winds. Some of the difficulties associated with 
TVD can be overcome by the use of active or passive aerodynamic surfaces, but these 
techniques encounter other difficulties especially with regard to size, location, control, 
and integration with TVD or reaction steering devices. It is apparent that important de- 
sign and performance problems are encountered when vehicle control is provided by 
either TVD, reaction, or aerodynamic means. Therefore, it is important to define real- 
istic steering requirements for solid-propellant rocket vehicles so that optimum control 
system configurations can be established. 

Minimum control requirements will be derived by assuming that the control force is 
obtained by thrust deflection. It should be noted, however, that the simplified analysis 
can be applied directly to active aerodynamic surfaces or reaction control devices. Be- 
cause these control devices are located ahead of the center of gravity, they tend to re- 
duce sideslip, which improves the accuracy of the simplified analysis. 

In the past, TVD requirements have been established by simulating the vehicle flight 
through a family of winds and noting the maximum thrust deflection angle used (ref. 1) . 
Most present-day vehicles use conventional control systems that utilize attitude and atti- 
tude rate feedback. They are designed to trim the vehicle during powered flight, that is, 
to maintain the nominal flight attitude even in the presence of winds . Until the develop- 
ment of solid-fueled engines, TVD angle was not a limiting factor and most control sys- 
tems were designed for trim deflection requirements. Solid-fueled engines pose prob- 
lems in obtaining TVD, which make it desirable to reduce the deflection requirements. 

The purpose of this report is to establish the TVD stability limits for a rigid-body 
configuration and a family of expected flight winds. A stable vehicle will be defined as 
one which can be returned to its desired flight path by appropriate control action when it 
deviates as a result of some disturbance. From the previous definition, the stability 
limit will be the minimum control action necessary to return the vehicle to its desired 
flight path after the disturbance subsides. The mathematical derivation of this stability 
limit is given in the analysis. Various assumptions are made to simplify the mathemati- 
cal model in order to give a closed-form solution. These assumptions and their effect on 
the results are discussed in the analysis. 

Validity of the theoretical results is established by detailed six degrees of freedom 
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(rigid body) computer simulations for two representative vehicle configurations: an in- 
termediate solid booster (260 in. (660 cm) solid) and a large solid vehicle capable of 
putting one million pounds of payload into 100 -nautical -mile (184-km) circular orbit. The 
required propulsion, aerodynamic, and weight data for these vehicles are based on a pre- 
liminary design study conducted at Lewis Research Center. The results include curves 
of theoretical and actual (detailed computer simulation) stability limits, and attitude er- 
ror convergence and divergence rates for TVD capability above and below the stability 
limit, respectively. The computer simulation and theoretical results are based on sim- 
plified wind shapes to reduce mathematical complexity, although the theory is applicable 
to any wind profile. 


ANALYSIS AND ASSUMPTIONS 

The equations of motion will be derived for a rigid-body configuration in order to 
simplify the analysis and to make the results easily applicable to vehicles of different 
size and aerodynamic behavior. The vehicle is assumed to have axial symmetry with 
negligible coupling between pitch and yaw planes. This allows independent pitch or yaw 
plane analysis. Thus the results of this report may not apply to vehicles with large wing 
sections or to vehicles having large coupling between pitch and yaw planes. 

The equations of motion are derived in a local horizontal, Cartesian coordinate sys- 
tem whose origin is located at the vehicle’s center of gravity. Aerodynamic moments 
are cancelled by thrust deflection. Because the duration of maximum flight winds is 
generally shorter than the atmospheric flight time, the analysis will be restricted to a 
time interval in which trajectory and vehicle parameters remain constant. Thus, the 
equations may be used to analyze vehicle response at any desired flight time. The re- 
sults in this report are given for maximum aerodynamic loading, which, in general, 
occur near maximum dynamic pressure. In order to obtain conservative control require- 
ments, maximum wind velocity and wind shear (the rate of change of wind speed with alti- 
tude) are assumed to coincide with maximum dynamic pressure. 

With these simplifying assumptions, there are two independent variables left in the 
problem, namely, wind velocity profile and control action against time. For statistically 
predictable winds (if the flight wind can be predicted in advance) , the maximum control 
requirements could be reduced by using more control action in the low-wind velocity re- 
gion than that required to trim the vehicle. This causes the vehicle to turn into the di- 
rection of high wind, thus reducing the angle of attack and, consequently, the control re- 
quirements in the high-wind region. In effect, the control requirement is averaged over 
a longer period of time, and the peak requirement is reduced. This control scheme has 
the possible shortcoming of building up too high an angle of attack due to excess control 
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before the high-wind-velocity region is encountered, wnich will limit the maximum control 
that can be used. With this consideration in mind, it is desired to determine how much 
can be gained by using maximum -control capability as soon as possible in the low-wind 
region, in other words, what percent of trim control is needed to maintain stability. The 
main purpose of this report is to obtain the theoretical minimum-control requirements 
applicable for most vehicle configurations. To this end, a general closed form stability 
boundary condition will be established. This has the form of an integral relation between 
the wind distribution and the necessary control action which will ensure a stable flight. 
The control system will be assumed to operate open loop; no attempt is made to design a 
closed-loop system. To illustrate how this generalized stability equation may be applied, 
the stability boundary will be derived for two different wind profiles and control actions. 


Derivation of General Stability Condition 

The following discussion includes the derivation of general stability boundary limit. 
To illustrate the usefulness of this result, specific boundary limits will be derived for 
rectangular and triangular wind profiles. The analysis is carried out for a rigid-body 
configuration with the assumption of uncoupled control planes; that is, the pitch and yaw 
planes do not interact. Figure 1 shows the basic vehicle configuration and the parameters 
employed in the analysis. The variables used are defined in appendix A. The vehicle is 
assumed to be aerodynamically unstable with thrust-vector control for stabilization. The 
reference trajectory is a zero angle -of -attack flight. The basic vehicle equations of mo- 
tion are given in appendix B. The principal analytic results can be derived from the mo- 
ment equation 


£ M CG= I e=N( a a-T7 c 5 (B2c) 

This equation applies for small deflection angles and a normal force that is proportional 
to the angle of attack. 

Nominal trajectory . - The nominal flight is a no wind, zero angle of attack (or zero 
lift) trajectory, which defines the nominal flight profile, and pitch attitude (© n ). The 
nominal flight parameters are 


“n-° 


a 


wn 


= 0 


4 



I 



Then, from equation (Blf), the nominal flight path angle is given by 


= © 


n 


n 


From equation (1) 


8 


n 



( 1 ) 


This nominal gimbal angle requirement is largest during pitch-over, that is, the time in- 
terval immediately after the vertical rise maneuver, which occurs before the high -wind 
velocity region. Because the gimbal requirements due to winds at maximum aerodynamic 
load conditions are much larger than the initial pitch-over gimbal requirement, this por- 
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tion of the flight poses no control problem. Also, the nominal control requirements in the 
high -wind region are very small, thus their contribution to the overall TVD angle is small. 
Therefore, the portion of the powered flight that can encounter stability problems due to 
deflection limitations is the high-wind velocity or maximum dynamic pressure region. 

Trim case. - The flight is trimmed if the aerodynamic moments are instantaneously 
cancelled by thrust deflection (this assumes an ideal autopilot). This will force the ve- 
hicle to maintain the nominal flight attitude (0 n ) even in the presence of winds: O t = © n - 
However, the vehicle position in an inertial coordinate system will deviate due to side- 
slip. Sideslip is the integrated effect of acceleration due to aerodynamic forces and the 
thrust component necessary to cancel the aerodynamic moments. From equation (Blc) 





Substituting equation (1) into this equation gives 


5 T- 6 n _ N *a = p 


On 


T l. 


( 2 ) 


The parameter B determines gimbal requirements for trim . One of the purposes of this 
report is, then, to determine how much this trim gimbal requirement can be reduced 
while still maintaining stability. The trim angle of attack using equation (Blf) is given by 

«r - - (r T - r n ) - % (3) 

If the variation in flight path angle is assumed to be negligible, the trim deflection is 
proportional to the angle of attack due to wind, thus equation (2) gives a simple formula 
to calculate trim deflection. 

Stability criteria . - In order to maintain a vehicle on some desired flight path, a 
certain amount of control action is required to counterbalance aerodynamic disturbances. 
A trimmed vehicle is by definition stable. The question that arises is what happens to a 
vehicle with a given control capability (5^) which is smaller than the trim control re- 
quirement for a given wind disturbance? There are two possibilities: The vehicle may 
lose control; that is, the attitude error (0 - © n ) increases indefinitely. Or the vehicle 
may return to the desired flight path after the disturbance ceases; that is, the attitude 
error increases, then goes to zero. In this report, the vehicle will be considered stable 
if it maintains the desired flight path or if it returns to the desired flight path after it has 
deviated due to some disturbance. The stability boundary is the minimum control action 
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necessary for stability divided by Ba . In general, the maximum trim deflection 

W y 

requirement is given by Ba™ . However, if the trajectory does not drift, then 

X y 1113^4 

Urj, = -a (as shown by eq. (3)), and the trim deflection requirement becomes 

1 y IIiaA W y lUaX 

Ba . The stability boundary obviously is a function of the vehicle parameters and 

w y mdx 

the wind duration and shape. The mathematical definition of stability is given by 

(© - ©J = 0 (4a) 

t-^oo 14 

From this equation the stability boundary can be obtained. 

The solution of (0 - © n ) contains exponential terms, whose coefficients are equal to 
the residues of the Laplace transform of (0 - 0 n ). In order for (0 - 0 R ) to have a stable 
solution, the residues of the unstable modes of (0 - 0 fl ) must all be zero; that is, 

lim (s - p.)#(e - 0 J = 0 for i = 1, 2, . . . , N (4b) 

s^ i 1 n 

where 0. are the unstable modes. This equation gives the general stability boundary; 
its application will be given in the following paragraphs. 

It can easily be shown that, if the - 0 n ) has the form 

•*(© - e n ) = ^[n x (s) JT(fl - 5 n ) + N 2 (s)jS?(a w )| 

where Nj, Ng, and D are polynomials in s, and /3. are the unstable roots of D, the 
stability condition can be expressed as 

- S ?(6 - 6 t ) = 0 

at s = i3. for i = 1, 2, . . . , N. This equation may be used if the trim deflection pro- 
file and the 0^ are available. 


Simplified Stability Analysis 

To simplify the moment balance equation (eq. (B3)), it is assumed that the flight- 
path angle is insensitive to winds and control action; that is, y = y n = y, p. With this 
assumption, the moment balance equation of appendix B becomes 
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( 5 ) 


s - p 


The finiteness of 6-6 and a assure that their transforms cannot have poles in the 

H W 

right half of the plane, thus, the only possible unstable mode is located at s = 0. Apply- 
ing the stability boundary condition of equation (4b) gives 


lim 

s— /3 


^(6 - 6 n ) = -B 


lim 

s H3 


^(« w ) 


(6a) 


or 



e -/3t (6 - 6 n ) dt = -B 




dt 


(6b) 


This equation is one of the important analytical results of this report. It gives an 
integral relation between any wind history and the corresponding gimbal-angle require- 
ment for stable vehicle control. To show how this equation is applied, two particular 
wind profiles will be considered with two different control schemes. 

Rectangular wind . - The wind is assumed to be a rectangular pulse of duration t 
and amplitude a . For control, assume that the thrust is instantaneously deflected, 

as soon as the wind is observed, to its maximum value 6 max . Sketch (a) gives a pic- 
torial representation of both wind and control. 
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Using equation (6b) yields 


6 


max 



e“ 0* dt = 


-Bff 


w, max 


f 

J 0 




dt 


or 

6 

_ max = l _ 6 -t/3 
Ba 

w, max 

The quantity -6 max /Bo' w max is defined to be the boundary limit. It is approxi- 
mately equal to the thrust deflection required for stability divided by the maximum trim 
deflection required. The formula is exact, as shown in the following equation, if the 
flight-path angle is unchanged from the nominal to the perturbed flight. 


yBR = 


1 - e” T ^ 


( 7 ) 


Therefore, 


y BR “ 


max 


Ba 


w, max 


Using equations (2) and (3), the boundary can be written as 


6 - 5 . 


y BR = 


n 


6 _ - 6 
T n 


( 8 ) 


The results of equation (8) are plotted in figure 2. For the vehicles considered, 0 is be- 
tween 0.4 and 0.6 radian per second. For large winds, the wind duration is at least 
10 seconds; therefore, t| 3 > 4. For these values of Tj3, the decrease in gimbal require- 
ments from the trim value is negligible or zero so that full trim control would be required 
to maintain stability when rectangular wind profiles are encountered. Fortunately, large 
rectangular wind profiles are not encountered. 




Figure 2 . - Theoretical stability boundary without flight- 
path angle effect. 


In order to obtain the attitude response to rectangular winds, the inverse Laplace 
transform of equation (5) must be evaluated. Define 


u(t - x) = • 


for 0 < t < x 
for x < t 


> x positive 


then 


6-6 =6 u(t) 
n max ' ' 


a w = a w,max[ u W - ^ - T )1 

Taking the Laplace transform and substituting into equation (5) gives 


p 


if(0 - 0 n ) = 


s 2 - P 


* 


±|6 + Ba (1 

2 s | max w, max' 


e -TS ) 


The inverse transform is given by 


10 



( 9 ) 


■ (1 - yBR> [ cosh t K;) ■ l ] u(t) • [ cosh T Kr ' ') ' l ] u (; • *) 


Triangular wind . - In this section, the wind is assumed to be triangular, which 
closely approximates synthetic winds often used in preliminary design studies. For con- 
trol action, assume that the vehicle is trimmed until time (t = a), then apply maximum 
control capability. 

At time t = a, the thrust deflection is equal to K times the maximum gimbal capa- 
bility. This corresponds to a dead-band in gimbal deflection of K6„„ . within which the 
vehicle is to be trimmed. Outside this band, maximum control is used, thus a is given 
by 


a 


max 


Ba 


w, max 



( 10 ) 


where 6 and a are the maximum deflection requirement and maximum wind 

IIlaA W y IllaX 

angle of attack, respectively. 

Referring to equation (7) , a becomes 


a 



( 11 ) 


Sketch (b) gives a pictorial representation of the control action and wind profile. 
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Using equations (2) and (3), the control action is given by 


or 


6-6 


n 


Ba T <» -Bff w 
6 

^ max 


for 0 < t < a 
for a < t 


6 - 6 n = -B[u(t) - u(t - a)]a w + 6 max u(t - a) 


where a is given by 
w 


a 


a = 
w 


w , max 
r 
2 


u(t)t - 2u ^t - — ^ |t - — ^ + u(t - r)(t - t) 


Taking the Laplace transforms and applying the boundary condition given by equation (6a) 
result in the stability boundary 



e -(r0/2)y BR K 


= e~ Tl3 (2e T/3 ^ 2 - l) 


( 12 ) 


The results of this equation are also plotted in figure 2. For r/3 = 4 and 100-percent de- 
flection dead band, the reduction in deflection requirement from the trim value is approxi- 
mately 31 percent. For smaller values of deflection dead band, the reduction is much 
greater. For a 0-percent deflection dead band, the maximum reduction is 62 percent of 
trim value. 

Substituting if (6 - 6 ) and Sf(ce) into equation (5) and taking the inverse Laplace 
n w 

transform give the attitude response to triangular winds in terms of the boundary condi- 
tion; that is, 
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Attitude error per angle of attack 


3- 0 1 


2.0 


l-Ol 


0 


Maximum deflection angle 

20 Percent above boundary limit 
10 Percent above boundary limit 
Boundary limit 

10 Percent below boundary limit 
20 Percent below boundary limit 




(c) Tp - 6. 

Figure 3. - Simplified analytic attitude error convergence and divergence rates. 
No flight-path angle effect; deflection dead band, 100-percent. 
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Attitude error per angle of attack 



(e) t|3 ■ 14. 

Figure 3. - Concluded. 



Attitude error per angle of attack 


1.0 


0 


Maximum deflection angle 

20 Percent above boundary limit 

10 Percent above boundary limit 

Boundary limit 

10 Percent below boundary limit 

20 Percent below boundary limit 




(a) ip = 2. 



(c) tp = 6. 

Figure 4. - Simplified analytic attitude error convergence and divergence rates. 
No flight-path angle effect; deflection dead band, 50 percent. 
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Maximum deflection angle 

20 Percent above boundary limit 

10 Percent above boundary limit 

Boundary limit 

10 Percent below boundary limit 

20 Percent below boundary limit 




(d) t|5 = 10. 


I 

t 

l 

i 

/ 



.60 . 80 1.00 1.20 1. 4( 

Time per wind duration, t It 


(e) xp ■ 14. 

Figure 4. - Concluded. 



Equation (13) is plotted in figures 3 and 4. 

Optimum control . - It is shown in appendix C by the use of calculus of variations that, 
if the wind time history is given, the gimbal angle against time history, which gives the 
smallest peak absolute value of gimbal angle, is the one where 6 - 6 n is constant. For a 
wind starting at time zero and control action applied at t = tj, using equation (6b), the 
optimum control is given by 

The factor e is the penalty for initiating a control action late, or the advantage that 
can be gained by starting early (t^ < 0) . The preceding equation gives the absolute mini- 
mum deflection that will maintain a stable flight. 

General Stability Analysis 

Appendix B develops the linearized differential equation for the flight-path angle. 

The resulting equations (B3) and (B6) form a complete set for the system. Combining the 
two equations results in equation (B8); that is, 


se(e - e n ) = [(s + bj) sr{6 - 6 n ) + b(s + b 2 ) *(og] (14) 

s 2 (s 2 - p 2 ) 

Applying the boundary condition given by equation (4b) to equation (14) gives another im- 
portant result of the analysis; that is, 

Z * (6 ■ , » ) = _BKi Z y(a '» ) (i5a) 


where 


0 + b 2 

K, = < 1 

1 0+bi 


(15b) 


Rectangular wind . - The wind profile and control action are assumed to be the same 
as in the previous section on rectangular wind. Using equation (15a), the boundary condi- 
tion becomes 
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^BR = K l (1 - e_T ^ 


(16) 


This equation shows that the deflection requirement is reduced if changes in the flight path 
angle are included. This was expected because it represents the amount of sideslip 
caused by the disturbance, which tends to reduce the angle of attack and consequently the 
deflection requirement. The amount of reduction depends on the particular vehicle con- 
figuration, because /3, b^, and bg are functions of the configuration. Generally speak- 
ing the reduction in deflection requirement would not be enough to preclude using full trim 
control just as in the case of no y effect, since Kj ^ 1. Attitude error convergence and 
divergence rates can be obtained in the same way as in the previous sections, by taking 
the inverse Laplace transform of equation (14) . 

Triangular wind. - The wind profile is assumed to be the same as in the previous 
section on triangular winds. Control action will also be the same; however, in this case, 

6 - 6 n will include the flight -path angle change from the nominal. The control action is 
given by 


6 - 6 n 


6 T- 6 n 


max 


for 0 < t < a 
for a < t 


From equation (14) 


Wt - 6 n ) - -B 


s + b„ 


s + b. 


-^(« w ) 


In order to obtain the if (6 - 6 n ), the inverse Laplace transform of the equation must 


first be evaluated. Then taking the Laplace transform of 


6 - 6 n - [u(t) - u(t - a)](6 T - 6 n ) + 6 max u(t - a) 


gives 


- S n ) = 


Bor 


w, max 

T 

2 


S + bq 


s^(s + bj) 


b,a 


1 r .. „ 1 . b 2 ' b l / e ^ 


2 2 yBR ^ ” ^ _ + 


s- - S b? V s +b l 


b. b x a v 
+ — + 


2 s Si 


-as 


(17) 


18 



Applying the boundary condition given by equation (15a) and triangular wind profile gives 


1 


" ~ y BR^ " *0 + ( b 2 b l) 




= K 1 e~ T ^(2e T ^ 2 - 1 ) 


(18) 


The results of equation (18) are plotted in figure 5 for 100- and 50-percent deflection 
dead bands. Results are presented for the two vehicle configurations discussed earlier 



Figure 5. - Theoretical stability boundary with and without flight-path angle effect for triangular wind. 


and compared with the case without flight -path angle effect given in figure 5. As it was 
expected, the boundary limit is reduced by including the flight -path angle change. For the 
case of 100 -percent deflection dead band and a r/3 = 4, the reduction is about 4 percent 
compared with the boundary with no flight -path angle included. The reduction is signifi- 
cant for larger values of tj3. 

The attitude response has been evaluated by taking the inverse transform of equa- 
tion (14) with the appropriate deflection and wind profiles. The resultant equation is 
given in appendix B (eq. (B9)) and is plotted in figures 6 to 9. These figures are dis- 
cussed in more detail in the RESULTS section. 

Optimum control . - It is shown in appendix C that, for a given wind profile, the gim- 
bal angle against time profile, which gives the smallest possible peak gimbal angle, is 
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the one in which 6 - 5 n = ±6 max - More specifically the optimal control function con- 
tains (N - 1) segments for which 5 = 8 or 5 = -6 where N is the number of 
real pole locations of ^0 - © n in the right half of the plane. It is also shown that, if one 
of the poles is much larger than the others (this is the case for the present problems) , 
the minimum value of 6 max is nearly the same as for the single-pole case, except for 
triangular winds with large values of r. 



Figure 6. - Analytic and actual attitude error convergence and divergence rates for vehicle 
configuration I. Flight-path angle effect included; deflection dead band, 100 percent. 
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(d) t(3 - 10. 

Figure 6. - Concluded. 
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Attitude error per angle of attack 



(c) t{3 = 6. 

Figure 7. - Analytic and actual attitude error convergence and divergence rates for vehicle 
configuration I. Flight-path angle effect included; deflection dead band, 50 percent. 


Attitude error per angle of attack 



(e) t|3 " 14. 


Figure 7. - Concluded. 



Attitude error per angle of attack 




Time per wind duration, t/i 


(c) x{3 = 10. 

Figure 8. - Analytic and actual attitude error convergence and divergence rates for 
vehicle configuration II. Flight-path angle effect included; deflection dead band, 
100 percent. 


24 



Maximum deflection angle 
10 Percent above boundary Analytic results 


— 10 Percent below boundary limit 

— « P ercen ‘ ab ° ve ]l m !H Actual computer run 

— 10 Percent below boundary limit J 
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(a) t(3 - 2. 
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(b) t| 3 a 4. 



.25 . 50 . 75 1.00 1.25 1.50 1.75 

Time per wind duration, t/x 


(c) Tp = 6. 

Figure 9. - Analytic and actual attitude error convergence and divergence rates for 
vehicle configuration II. Flight-path angle effect included; deflection dead band, 
50 percent. 




COMPUTER SIMULATION 


To show the accuracy of previously obtained equations, a six degrees of freedom tra- 
jectory program was used to obtain actual stability boundaries. Computer simulations 
with detailed aerodynamic, propulsion, and inertia calculations were made for two vehicle 
configurations, the 260-inch (660-cm) solid booster (configuration I) and a large solid ve- 
hicle (configuration n) capable of delivering 1 million pounds (453 600 kg) of payload to 
100-nautical-mile (184-km) orbit. Weight, inertia, and aerodynamic force coefficients 
have been generated at Lewis. Both vehicles were targeted to a 100-nautical-mile circu- 
lar orbit with the constraint that the dynamic pressure will not exceed 950 pounds per 
square foot. The booster portion of flight, or aerodynamic flight, is constrained to a 
zero lift path to establish the nominal flight profile 0 , which was curve fitted for the 
nonnominal trajectory simulations. Pertinent trajectory and vehicle parameters used to 
obtain theoretical boundary values are given in table I for both configurations at maximum 


TABLE I. - NOMINAL TRAJECTORY PARAMETERS AT MAXIMUM 
DYNAMIC PRESSURE 



Configuration 


I 

n 

Weight, W, N 

1.48X10 6 

1.393X10 7 

Moment of inertia, I, (kg)(m ) 

0. 225x10® 

1.033X10 9 

Normal force per angle of attack, N, N/rad 

0.667X10 6 

0.562X10 7 

Total thrust, T, N 

2. 495X10 6 

3.334X10 7 

Aerodynamic moment arm, Z , m 

9.45 

28.78 

Control moment arm, l , m 

23.56 

50.27 

Axial force, F^, N 

0.807X10 5 

0.962X10 6 

Nominal relative velocity, V , m/sec 

538.81 

614.23 

Nominal flight attitude, © deg 

44.83 

37.52 

Time, t, sec 

66 

85 


dynamic pressure conditions. Wind profiles correspond to those used in the analysis. 
Synthetic winds can be closely approximated by a triangle of appropriate time duration. 
Real winds are more difficult to approximate by rectangular or triangular shapes. How- 
ever, reasonable results can be obtained by approximating the real wind with a triangle 
of the same area and maximum velocity. 

The control system was programmed to apply trim control until the deflection re- 
quired reached K percent of maximum, then use maximum allowable control. 
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RESULTS 


Boundary values given by equations (8), (12), and (17) are plotted in figures 2 and 5. 
In figure 2, the percent of trim deflection required for stability is shown for rectangular 
and triangular wind profiles with no flight -path angle effects included. Boundary values 
for 100-, 75-, 50-, 25-, and 0-percent deflection dead bands are presented for triangular 
winds. As is seen from this figure, for rectangular wind profiles, the vehicle must be 
trimmed at all times, because 98 percent of trim deflection is needed to maintain sta- 
bility for a typical value of t/3 = 4. For the same value of t/3 in the case of triangular 
winds, only 69 percent of trim control requirement is needed for a 100 -percent deflection 
dead band. This value can be further reduced by anticipating the wind to behave in a cer- 
tain way. Obviously, if the wind velocity is assumed to continue to increase, the best 
control action would be to use maximum control capability as soon as the wind is ob- 
served. In this way, the control requirement can be reduced to 38 percent of trim value 
for zero-percent deflection dead band. This procedure defines the minimum control re- 
quirement attainable. Unfortunately, it requires accurate wind prediction. 

A more realistic control procedure would be to trim the vehicle until the wind angle 
of attack exceeds some fixed limit, then apply maximum control. A 50-percent deflection 
dead band which is switching at half maximum control capability, reduces the deflection 
requirement to 56 percent of trim for t/3 = 4, or more for long wind duration. Figure 5 
gives a comparison between boundary limits with and without flight path angle effects for 
different values of dead band. As was expected, there is large reduction in control re- 
quirements due to flight -path angle effect for long wind duration because of sideslip or 
drift velocity, which tends to reduce the angle of attack. In figure 5(a), for r/3 = 4, the re- 
duction is about 4 percent, and it increases to 24 percent for tj3 = 20. Figure 5(b) is plot- 
ted for a deflection dead band of 50 percent, and the reduction in control requirement is 
negligible because the angle of attack is reduced and the sideslip is essentially eliminated. 

The results given in figure 5(b) indicate that approximately 56 percent of trim con- 
trol capability will maintain stability, although the vehicle will deviate from the nominal 
flight profile. The amount of attitude error, for triangular wind profiles and no flight- 
path angle effect included, may be obtained from figures 3 and 4, for different values of 
control dead band. The curves represent attitude error for control capability 10 and 20 
percent above and below the boundary value. For the 260-inch (660-cm) solid- fueled 
rocket with /3 = 0. 53 reciprocal second and a wind duration of 7. 5 seconds, will 
equal 4. Assuming 100-percent deflection dead band and with a control capability 10 per- 
cent above the boundary, the maximum attitude error is 0. 16° per wind angle of attack. 
For a 5° angle of attack, the vehicle attitude would deviate a maximum of 0.8° from the 
nominal. The attitude error goes to zero at t/r =* 1, or t =* 7. 5 seconds. The attitude 
error continues in the negative direction after passing through zero. For the real auto- 
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Maximum deflection per maximum trim deflection Maximum deflection per maximum trim deflection 
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Tp 

(a) 100-Percent deflection dead band. (b) 50-Percent deflection dead band. 


Figure 10. - Theoretical and actual stability boundary for vehicle configuration I. 
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(a) 100-Percent deflection dead band. (b) 50-Percent deflection dead band. 


Figure 11. - Theoretical and actual stability boundary for vehicle configuration II. 
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pilot, however, the engine deflection would be decreased to the value required to maintain 
control after the attitude error has gone to zero. Therefore, the negative portions of 
figures 3 and 4, and 6 to 9 should be neglected as far as attitude error is concerned. The 
shape of these curves, however, give valuable information about the rate of convergence 
to the nominal attitude. The curves show rapid convergence and divergence rates above 
and below boundary value conditions, respectively. For triangular winds, the attitude 
error in most cases goes to zero before the wind subsides with a control capability 
10 percent above boundary value . The exceptions are short wind gusts or rectangular 
wind distributions. In these cases, the time required to eliminate attitude errors is of 
the order of 1 to 2 times the wind duration. 

Figures 10 and 11 give a comparison of actual and theoretical boundary limits for the 
two vehicle configurations. In both cases, there is good agreement with the actual six 
degrees of freedom computer simulations. In the analysis, it was assumed that the tra- 
jectory parameters will remain constant. However, in the actual flight they do vary, and 
their variations limit the time interval to which the theoretical results may be applied. 
This limitation could be partially circumvented by doing a point by point calculation along 
the trajectory with the appropriate wind for each segment. However, in general, it is 
sufficient to obtain deflection requirements at maximum aerodynamic load conditions. 

A comparison of actual computer results and the theoretical results with flight-path 
angle included are given in figures 6 to 9. These figures show that the theoretical re- 
sults are conservative; that is, the convergence rate is slower while the divergence rate 
is faster than the actual computer results. 


CONCLUSIONS 

The equations presented and the analytic results derived give a concise and complete 
description of an open-loop minimal control requirement to maintain stability for a gen- 
eral rigid-body vehicle configuration. The formulas are easily applicable, with the added 
advantage that the analysis can be carried out at any particular flight time. The results 
are given for the pitch plane, but, with minor modifications, they can also be applied to 
the yaw plane. 

General stability boundary conditions are presented which do not depend on the ve- 
hicles considered in this report, and can be easily applied to any vehicle configuration by 
computing a single vehicle dependent parameter. It is also shown (for the vehicle config- 
urations considered) that these boundary conditions give slightly conservative results. 

For more accurate results, several additional vehicle- and trajectory-dependent param- 
eters must be calculated and a transcendental equation must be solved. 
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The results obtained show that, with good control-system design, the deflection re- 
quirements can be reduced to about 56 percent of trim deflection. At the same time, be- 
cause of the method of control, the maximum angle of attack and bending moments can 
also be substantially reduced. 

Rapid convergence and divergence rates are also shown for control capability above 
and below the given stability boundary, respectively. This implies rapid convergence to 
the nominal flight profile after deviation from it due to some disturbance, if the control 
capability is only slightly above the stability boundary value. Although the convergence 
rate is rapid, in designing a closed-loop control system one should not get too close to 
the boundary limit because, if the system becomes unstable for any reason, the diver- 
gence is also rapid. 

A comparison of actual (detailed computer simulation) and theoretical attitude errors 
indicates that the theory gives somewhat conservative results; that is, the theoretical 
convergence rate is slower while the divergence rate is faster than the actual. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 31, 1968, 

125-17-05-01-22. 
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APPENDIX A 
SYMBOLS 


a 

control switching time defined 

t 

time, sec 


by eq. (10), sec 

*8 

switching time defined in 

B 

constant defined by eq. (2) 

appendix C, sec 

b V b 2 

constants defined in appen- 
1 

u(t) 

unit step function 

C. 

1 

dix B, sec 

constants defined in appen- 
dix B, sec 1 

normal force coefficient per 
angle of attack, rad - '*' 

V 

relative velocity including 
wind, m/sec 

dc n /do- 

V R 

V w 

relative velocity, m/sec 
wind velocity, m/sec 

D i 

defined in eq. (Cl) 

W 

gross weight, N 

E 

Wierstrass E variable 

y BR 

boundary limit 

f a 

axial force, N 

a 

angle of attack, rad 

f i 

constraint equations defined 

a w 

wind angle of attack, rad 


in appendix C 

0 

vehicle response parameter, 

g 

gravitational acceleration, 


sec -1 

I 

m/sec 

moment of inertia, (kg)(m^) 

*i 

unstable roots, sec - * 

J 

function to be minimized de- 
fined in appendix C 

y 

5 

flight path angle, rad 
deflection angle, rad 

K 

multiplier (percent deflection 

0 

flight attitude, rad 


dead band) 

X. 

1 

undetermined Lagrange multi 

K 1 

constant defined by eq. (16) 


pliers 

Z a 

aerodynamic moment arm, m 

P 

Q 

air density, kg/m 

Z c 

control moment arm, m 

a 

real variable defined in ap- 

M cg 

total moment about center of 


pendix C 

gravity, N 

T 

wind duration, sec 

m 

N 

gross mass, kg 

normal force per angle of 
attack, N/rad 

Subscripts: 

max maximum 

S 

2 

vehicle base area, m 

n 

nominal 

s 

Laplace variable, sec”^ 

T 

trim 


T total thrust, N 


31 



APPENDIX B 


EQUATIONS OF MOTION 

The basic vehicle equations of motion will be derived for a rigid-body configuration. 
Only motion in the pitch plane is considered. The equations of motion and the vehicle’s 
control system are treated as uncoupled, insofar as the vehicle's three axes of rotation 
are concerned. The equations are given in the pitch plane, but the yaw plane analysis can 
be done in a similar fashion. In order to apply the following equations to the yaw plane, 
the only changes that have to be made are to set © n = 90° and g cos © n /V n = 0 in equa- 
tion (4). The basic vehicle configuration is given in figure 1, and the reader is referred 
to appendix A for the definition of variables and symbols used. 

Basic equations of motion: 

^ F|| y = mV = T cos (0 - y + 6) - F A cos (0 - y) - No- sin (0 - y) - W sin y (Bla) 

I F lV = mVy = T sin(© - y + 5) - F^ sin(0 - y) + No- cos (0 - y) - W cos y (Bib) 

y^M CG = I© = NZ a a - Tl c sin 5 (Blc) 

1 9 dc 

N = — pVp — ^ S (Bid) 

2 R da 

V w 

sin » w = — — sin y (Ble) 

V R 

a = © . y _ a w (Blf) 

Using small-angle approximations for a, a and 6, the equations become 

^Tf h y = mV = T [cos (0 - y) - 5 sin (© - y)] - F A cos (0 - y) - Ntt sin (0 - y) - W sin y 

(B2a) 

Z F 1V = mVy = T[sin (0 - y) + 5 cos (© - y)J - F A sin(© - y) +Na cos (© - y) - W cos y 

(B2b) 
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Z M Co'® sNi a'” T y 


w 

a = sin y 

w v 

V R 


a = © - y - a 


w 


(B2c) 

(B2d) 

(B2e) 


To derive the principal analytic results of this report, it is necessary to express the 
attitude error in terms of the wind and control action. The attitude error is given by 
0 - © n or, from equation (B2c), 

I(e - e n ) = N; a (o. - % ) -T* c (6-5 n ) 


Substituting for (a - ar ) from equation (B2e) into the preceding equation and taking the 
Laplace transform give 


- © n ) = 


|^(r-r„) + |[i?(8-5n) + B ^(«w3} 


(B3) 


where 


/3 


2 




(B4) 


In order to eliminate jS f(y - y ) from equation (B3), equations (Bla) and (Bib) must be 
linearized about some nominal operating point. 

Linearized equations for the flight-path angle . - In the following analysis, the same 
symbols will be used for the linearized variables as for the actual variables. Keeping 
the same symbols simplifies the notations and makes the equations less confusing. Lin- 
earizing equations (B2a) and (B2b) gives 

mV = -T sin (© - y + 6) (0 - y + 6) + F. sin (© - y)_(0 - y) - Nor cos (0 - y)„(0 - y) 

n n n n 

- W cos y n y - Na sin (0 - y) n 
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m(Vy n + V n y) = T cos (© - y + 6) n (0 - y + 5) - cos (0 - y) n (0 - y) 


- Na n sin (© - y) n (0 - y) + Nff cos (0 - y) n + W sin y n y 


but 6 n « 0 n , o? n = 0, and 0 R = y R ; therefore, 


mV = -W cos 0 v 
n' 


m(Vy + V y) = T(0 - y + 5) - F. (0 - y) + Na + W sin 0 n y 
n n a 


Substituting for V n and y n from equations (Bla) and (Bib) and taking the Laplace 
transform give 


y(y) = 


s 


mV n 2 T - F A - W sin 0 n /g cos 0\ 2 


if |(T - F a )0 + T5 + Na] 


mV. 


n 


n 


Substituting for a from equation (B2e) gives 

sif|(T - F a + N)0 + T5 - NckJ 


nr) = 


mV 


n 2 T-F a+ N-W sin © n 


^g cos © 


s + 


s + 


mV 


n 


(B5) 


Because 


(g cos 0 \ T - F a + N - W sin 0 
2. « A n (s| 


n 


mV. 


n 


there are two real roots, one of which is located very close to the origin. To simplify 
the analysis, the smaller root is assumed to be at the origin. Then 

Se(y) = — — if f(T - F a + N)0 + T5 - No 1 (B6) 

mV s + C n L A WJ 

n l 
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where 


C l = 


T - F. + N - W sin ©„ 
A n 


mV, 


n 


Substituting equation (B6) into equation (B3) gives 




&{© - © ) = - J*_ |(s + b 1 )^f(6 - 6) + B(s + b 
P(s) L n 


where 


P(s) = s 3 + CjS 2 - p 2 s + C 2 


g sin © n 
C 2 = P V 


n 


b i = 


mV 


n 


T - F A + N 4 + - w sin € 


b 2 = — — (t - F. - W sin ©J 
4 mV n 


n 


Since C 1 and C 2 are much smaller than /3 2 , the polynomial 


P(s) = (s - 


/ 3 2 /L 


s 2 + [ c i+ -2 s -^ 2 




Also, |s| [ Cl + (C 2 /^ 2 )]« /3 2 and the two roots are located at 

is very small [(C 2 / fi 2 ) « O.Olj and can be assumed to be 
P(s) becomes 


zero. 


P(s) » s(s 2 - j3 2 ) 


'2>^w>] 


(B7) 


P(s) can be factored as 


s « ±/3. The root C 2 //3 2 
With these assumptions 
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therefore, 


S?(Q - 0 n ) ? [(s + b 1 )^f’(6 - 6 n ) + B(s + b 2 ) ( B8 ) 

s(s 2 - 13 2 ) 

The inverse Laplace transform of the attitude error response to triangular wind and 
trim deflection profile with dead band is given as follows: 


0 - 0 , 


n 


-a 


w, max 


_2 

n 3 


biT 


1 + a(b 2 - b,) - y br (1 - K) 


sinh r/8 (— - - t/3 (— - — 

\T r) \T T, 


+ — cosh t/3 (— - — \ 

fi \r 




'B(b 2 - bj) 


b; 


-b^a 

e 1 + abj - 1 


TP_ 

2 


Y br (1 - K) 


cosh T/S 


c- 



- 2 -<sinh t/3 
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APPENDIX C 


EQUATIONS FOR OPTIMUM CONTROL 

The problem of finding the optimum deflection profile for a predictable wind profile 
can be stated mathematically as follows. It is desired to find the function 6(t), with the 
smallest possible maximum value, that satisfies a set of definite integral constraints. 

The constraints are given by equation (15a) 

The problem is formulated mathematically as a Bolza problem of the calculus of 
variations; namely, it is desired to minimize 6 max > where 6 max is the maximum value 
of 6(t) subject to the following constraints: 

/ “-fit -jfi.t 

e 1 6(t)dt + BK 1 / e 1 a w (t)dt 

u •'0 


-f 


-At 

e 1 6(t)dt 


D i 


i = 1,2,3, 


N 


(Cl) 


5 is defined by the additional constraint 
max 


% + l - 62 - 6 Lx + ^ 


(C2) 


where a is a real variable. The constraint equations are now adjoined to the function to 
be minimized in the usual manner, thus the function to be minimized J is 

N+l 

j=ZVi* w 

i=l 


where the A.. are the undetermined La Grange multipliers. The Euler -La Grange equa- 
tions to determine the optimum solution are 


A /A*_\ = i!L i = 1, 2, . . . , N 
dt \dX i J 6X i 
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where 


X x = 6 
X 2 = a 


Then 


/ ; x i e + 2 a N+ 1 5 ~ 0 
i=l 


(C3a) 


2 X N + l a 


= 0 


(C3b) 


From equation (C3a) 


N 


5 = -- 


2 A n 


N+l i=l 




In order for 6 to be bounded, Ajj + ^ must not be zero. Therefore, from equation (C3b), 
a must be equal to zero. Using this result and equation (C2) gives 


5 


± ^max 


(C 4) 


The choice of sign in equation (C4) is determined by the Wierstrass E-test, which gives 


E = J(3T ) - J(X) + ^L-(X - X* ) s 0 

ax 


for a minimum. For this problem 


E =12 v ( 6 - 6 > 

i=l 

where 5 is the minimum value and 5* the allowable variation; in this case, 6* = -5. 
Then 
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A -fit 

2 > , X^ 1 6 < 0 
i=l 


Because the X., are constants, the preceding function can have at most N - 1 zeros. 

These zeros correspond to the switching points from 8 = 6 to 6 = -6 and 

max max 7 


6 = 5 


max 


N 




if 2 v 1 - 0 

i=l 


6 = -5 


max 


V' -ft 

if \ e 1 > 0 
i=l 


Example: Suppose there are two integral constraints. Then the switching function 
becomes 


-/3,t -j8„ 1 

X^e 1 s + Xge s 


which gives a switching time of 


= -A) 

f 2'fl \ X 2 / 

The solution is obtained by choosing 6 and t such that equation (Cl) is satisfied. 

nidx s 


S max l 1 - 2e "' 3l ‘ S ) =0l D l 




2 u s 


“maxi 1 " 26 “V*i D 2 


Combining these equations gives 


0 oD 


2 2 


max 


-0 2 t e 
1 - 2e * £ 


(C5a) 
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(C5b) 


2e 


-J3 


1*8 


2 ^1 D 1 e " /3 2 t s 
^2 D 2 


1 - 


^2 D 2 


Let 


x = 



2*8 


y = 


^i D i 

^ 2 D 2 


Then equation (C5b) becomes 


2x 1 - 2xy = 1 - y 


(C6) 


If x is approximately zero or one, then equation (C5a) shows that the stability boundary 
is nearly unchanged from the single-pole solution. For most launch vehicles, /3g is at 
least a factor of 10 larger than /3j. If this is the case, the solution to equation (C6) for 
0 < y < 1 is 


-M 2 " 1 


(C7) 


Therefore, for 0 ^ y ^ 1, x is never larger than 0.001. For y = 1, equation (C6) has 

two solutions: x = 0 or x = 1. The two solutions correspond to t = 0 or t = «> and 

s s 

6 max = (3 2 D 2’ and are > therefore, physically equivalent. For y » 1, the solution to 
equation (C6) has the form 


x-yjLi 


2y 
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The total solution of equation (C6) is shown in sketch (c) 



From the preceding discussion and sketch, the single-pole stability boundary solution 
can be expected to be valid up to y = 1 . The value of y must now be related to t for 
triangular winds, in order to predict the range of validity of the single-pole solution. 

For a triangular wind duration r, 


D i = 


1 - e 


-r^/r 


i = 1.2 




and 


y 




(C8) 


In the limit 
approaches 
r at which 


as t approaches zero, y approaches which is less than 0. 1. As t 

», y goes to or ^ eas ^ *0- Therefore, there exists some value of 

y = 1 . From equation (C8) , the solution can be expressed as 


-T/V2 

e 



(C9) 


If jSg = 0. 5 and = 0.01 (the approximate value for the vehicles considered herein), 
then t = 28 seconds for y = 1. Thus, the single-pole approximation can be expected to 
be valid up to r = 28 seconds (or rj3 = 14) for a triangular wind. For longer wind dura- 
tion, the value of 6 max may increase substantially. 
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